Scaling exponents for fracture surfaces in homogenous glass and glassy ceramics 
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We investigate the scaling properties of post-mortem fracture surfaces in silica glass and glassy 
ceramics. In both cases, the 2D height-height correlation function is found to obey Family- Viseck 
scaling properties, but with two sets of critical exponents, in particular a roughness exponent £ ~ 
0.75 in homogeneous glass and £ ~ 0.4 in glassy ceramics. The ranges of length-scales over which 
these two scalings are observed are shown to be below and above the size of process zone respectively. 
A model derived from Linear Elastic Fracture Mechanics (LEFM) in the quasistatic approximation 
succeeds to reproduce the scaling exponents observed in glassy ceramics. The critical exponents 
observed in homogeneous glass are conjectured to reflect damage screening occurring for length- 
scales below the size of the process zone. 

PACS numbers: 62.20.Mk, 46.50.+a, 68.35.Ct 



The morphology of fracture surfaces is a signature of 
the complex damage and fracture processes occurring 
at the microstructure scale that lead to the failure of a 
given heterogeneous material. Since the pioneering work 
of Mandelbrot , a large amount of studies have shown 
that crack surface roughening exhibits some universal 
scaling features: Fracture surfaces were found to be self- 
afhne over a wide range of length scales, characterized 
by a universal roughness exponent £ w 0.8, weakly de- 
pendent on the nature of the material and on the fail- 
ure mode (see e.g. for a review). Very recent stud- 
ies showed that a complete description of the scaling 
properties of fracture surfaces calls for the use of the 
two-dimensional (2D) height-height correlation function. 
This function was observed to exhibit anisotropic scaling 
properties similar to the Family- Viseck scaling £| pre- 
dicted in interface growth models characterized by 
three critical exponents independent to some extent of 
the considered material, the loading condition and the 
crack growth velocity. 

The origin of the scaling properties of fracture surfaces 
is still debated. Hansen and Schmittbuhl [(| suggested 
that the universal scaling properties of fracture surfaces 
are due to the fracture propagation being a damage co- 
alescence process described by a stress-weighted perco- 
lation phenomenon in a self-generated quadratic dam- 
age gradient. Bouchaud et al. 7] proposed to model the 
fracture surface as the trace left by a crack front mov- 
ing through randomly distributed microstructural obsta- 
cles - the dynamics of which is described through a phe- 
nomenological nonlinear Langevin equation, keeping only 
the terms allowed by the symmetry of the system. Finally, 
Ramanathan et al used Linear Elastic Fracture Mechan- 
ics (LEFM) to derive a linear nonlocal Langevin equation 
within both clastostatic 8] and clastodynamic 9] approx- 
imation. All these approaches succeed to reproduce scale 
invariant crack surface roughness in qualitative - but un- 
fortunately not quantitative - agreement with the exper- 
imental observations 0. 

The universality of the roughness exponent was found 



to suffer from several exceptions: Metallic surfaces inves- 
tigated at the nanometer scale were found to exhibit self- 
affine scaling properties, but with a roughness exponent 
significantly smaller than 0.8, closer to 0.4-0.5 [ToL Hl|. 
This was first interpreted as a kinetic effect similar to 
the one expected for a moving line close to its depin- 
ning transition 11] - the small (resp. large) scale rough- 
ness exponent 0.5 (resp. 0.8) corresponding to effective 
quenched noise (resp. thermal noise) |l2^ . The relevance 
of such an interpretation was later questioned since no 
small scale £ ~ 0.4 — 0.5 roughness exponent was ob- 
served for nano-resolved fracture surface of Silica glass 
broken under stress corrosion with crack growth velocity 
as small as the picometer per second j3j- Furthermore, 
recent experiments reported similar values £ ~ 0.4 — 0.5 
at large length-scales in sandstone |l3j | , artificial rock [bH 
and glassy ceramics [l5|. In this latter case, the rough- 
ness exponent was found to be independent of the bead 
size, the porosity, the transgranular/intergranular nature 
of the failure mode and the crack growth velocity. This 
suggests the existence of a second universality class for 
failure problems. 

The series of experiments reported here were designed 
to uncover the origin of these two distinct universality 
classes and focus more specifically on the range of length- 
scales over which the scaling properties are observed. 
Two materials are investigated: homogenous glass and 
glassy ceramics made of sintered 100 /im glass beads. 
In both cases, the fracture surfaces are found to exhibit 
Family- Viseck scaling properties but with two different 
sets of critical exponents, in particular £ ~ 0.75 for ho- 
mogenous glass and £ ~ 0.4 for glassy ceramics. The 
range of length-scales over which these two scalings are 
observed are shown to be below and above the size of 
process zone, respectively. Using LEFM, we show that 
the crack roughness development can been described as 
an elastic string with nonlocal interactions creeping in a 
2D random medium - the spatial coordinate along which 
the crack globally grows playing the role of time. This 
approach allows to account quantitatively for the value 



2 



of the observed critical exponents in the case of glassy 
ceramics. The role of damage in the case of homogeneous 
glass is finally discussed. 

Experiments. - In all the following, the reference frame 
(e x , e y ,e z ) is chosen so that e x , e y and e z are parallel to 
the propagation, loading and crack front directions re- 
spectively. Fracture surfaces in amorphous silica were ob- 
tained for various growth velocity ranging from v = 10~ n 
to v = 10~ 4 m/s using the procedure described in Ref. 
@, H3]. Their topography was then measured through 
Atomic Force Microscopy with in-plane and out-of-plane 
resolutions estimated to be 5 nm and 0.1 nm respectively. 
The resulting images are 1024 x 1024 pixels 2 and repre- 
sent a square field of 1 x 1 ^m 2 . The scaling properties of 
the fracture surfaces were analysed using the procedure 
discussed in Ref. [3j: First, the ID height- height correla- 
tion function Ah(Az) =< (h(z + Az, x) - h(z, x)) 2 > 1/2 
along the z direction and Ah(Ax) =< (h(z, x + Ax) — 
h(z, x)) 2 > x / 2 along the x direction were computed. Both 
Ah(Az) and Ah(Ax) were found to exhibit power-law 
behaviours, characterized by exponents C — 0.75 and 
(3 ~ 0.6 respectively, extending up to length scales £ 2 
and £ x respectively (Fig.^ Inset). Second, the 2D height- 
height correlation function was computed (Fig. ^| and 
was shown to follow a Family- Viseck scaling jj]: 
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FIG. 1: Inset: ID height-height correlation function calculated 
along the propagation direction z and the crack front direction 
a; on a fracture surface of silica glass. In this experiment, the 
crack growth velocity was v — 10 -10 m/s. The straight lines 
are power law fits (see text for details). The vertical dot-dash 
line sets the cutoff lengths £ 2 ~ 70 nm and £ x ~ 50 nm of the 
self-amne regime in the z and x directions respectively, (b): 
The insets show the 2D height-height correlation functions 
A/iAa;(Az) corresponding to different values of Ax vs Az. 
The data collapse was obtained using Eq. with exponents 
C ~ 0.75, (3 ~ 0.6, and z = £//3 ~ 1.2. 



Ah oc Ax/V(a^), /(«) 



if u < 1 
if u > 1 



(1) 



with C ^ 0.75, (3 ~ 0.6 and z = C//3 ^ 1-2, for length- 
scales Az < £ z and Ax < £ x , in agreement with previous 
studies reported in These exponents £, (3 and z were 
shown to correspond to the roughness, growth and dy- 
namic exponents defined in interface growth problems. 

Then, we investigated the influence of the crack growth 
velocity v on the scaling properties of the post-mortem 
fracture surface. The critical exponents (, (3 and z do 
not show any noticeable dependence on v. On the other 
hand, the cutoff length £ was observed to decrease slowly, 
as the logarithm of v (Fig. |2J). For the smallest value 
of v, ranging from 10~ n to 10~ 9 m/s, we were able to 
observe in real-time, at the nanometer scale, the crack 
propagation during the specimen failure |18|, At these 
scales, it was shown that the deformation fields does not 
fit with the linear elastic predictions over a fairly large 
region (~ 100 nanometers) at the crack tip [lti UtI fl8| . 
This zone is thereafter referred to as the damage zone 
or the process zone. The variation of its size R c with 
respect to the crack velocity v are presented in the in- 
set of Fig. |21 First, R c is found to be larger, but of the 
same order of magnitude, than £. Second, R c , like £, is 
observed to decrease as the logarithm of v. This leads 
us to conjecture that the damage zone size R c is the 
relevant length-scale that set the upper cutoff length £ 
that limits the scaling given by Eq.^with the exponents 
{C ^ 0.75,/? ~ 0.6, z = C//3 ~ 1.2}. At these length- 
scales, the material cannot be identified with a coarse- 
grained equivalent linear elastic medium, which explains 



the failure of existing models [8j derived from Linear Elas- 
tic Fracture Mechanics (LEFM) to reproduce the critical 
exponents {( ~ 0.75, (3 ~ 0.6, z = Q/(3 ~ 1.2} observed 
exp er iment ally. 
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FIG. 2: Variation of the cutoff length £ (measured along z) 
as a function of the crack growth velocity V . The axes are 
semilogarithmic. The straight line corresponds to a fit £ oc 
log(w). Inset: Variation of the size of the damage zone R c 
(measured along x) as a function of the crack growth velocity 
v (see text for details). The axes are semilogarithmic. The 
straight line correspond to a fit R c oc log(u) 

We now examine the scaling properties of fracture sur- 
faces in glassy ceramics made of sintered beads of Silicate 
glass with diameter d ranging between 104 — 128 /im (see 
Ref. Q for details). In this class of materials, the size 
of the process zone observed in the vicinity of the vari- 
ous (micro-)crack tips are expected to be of the order of 
100 nm as in homogeneous glass, while the microstruc- 
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ture scale is set by the mean bead diameter at a length- 
scale three order of magnitudes larger. Figure presents 
the 2D height-height correlation function. As for homo- 
geneous glass, this function is found to follow the Family- 
Viseck scaling |4J given by Eq. ^ but with a different set 
of critical exponents {£ ~ 0.4,/? ~ 0.5, z = £//3 ~ 0.8}. 
This set of exponents was found to be independent of the 
porosity, the bead diameter and the crack growth veloc- 
ity. 
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FIG. 3: The insets show the 2D height-height correlation func- 
tions Ah Ax ( Az) corresponding to different values of Ax vs Az 
for a fracture surface of glassy ceramic made of glass beads 
with a porosity $ = 6%. The data collapse was obtained using 
Eq.EJwith exponents £ ~ 0.4, (5 ~ 0.5, and z = C//3 ^ 0.8. 



Intepretation. - Since the scaling properties of glassy ce- 
ramics were observed at length-scales larger than the pro- 
cess zone size, it is natural to attempt to interpret these 
scaling properties within LEFM framework. We restrict 
the following analysis to the case where (i) the elastic 
properties can be considered as homogeneous and (ii) the 
crack speed is small compared to the speed of wave prop- 
agation so that quasi-static approximation is relevant. In 
an isotropic elastic material, near the crack front, the 
stress field can be written fill: 



^M(f) + K n m {6) + KujW^O)}, 



where r is the distance from the crack front, 8 is the angle 
with respect to the local direction of crack propagation, 
Sy are universal functions of 9, and K[, Kn and Km 
are the mode I (tension), mode II (shear) and mode III 
(tear) intensity factors respectively. In an ideal homoge- 
neous elastic material under tensile Kj loading, the crack 
front would remain straight and would propagate within 
a given plane. But heterogeneities of a disordered mate- 
rial like glassy ceramics induce both in-plane distortion 
f(z, t) and out-of-plane distortion h{z, x = f(z, t)) of the 
crack front. In turn, these deviations from straightness 
induce perturbations S Ki, SKu and SKm in the local 
loading of the crack front. To linear order, SKi depends 
only on / [2(j while SKjj and 5Km are functional of h 
only |21| . Since we are primarily interested in the fracture 
surfaces, we focus on the interrelations between 5Ku and 



SKju and h, returning in the end of the paper to a brief 
discussion of the in-plane crack roughness. 

The path chosen by a crack propagating in an elastic 
isotropic material is the one for which the local stress field 
at the tip is of mode I type (" criterion of local symme- 
try") [illHIll In other words, the net mode II stress 
intensity factor should vanish in each location z along 
the crack front. Three contributions should be taken into 
account in the evaluation of 5Ku. The first contribution 
is due to both inevitable imperfections in the loading 
system or in crack alignment and the heterogeneous na- 
ture of the material. This contribution is modelled by a 
quenched uncorrelated random field SK H = Kji](z, x, h) 
written, for sake of simplicity, as the sum of two uncorre- 
lated random fields SK H — K®{r) q (z, h) + rjt(z, x)}. The 

(2) 

second contribution 8K U arises from the coupling of the 
singular mode I component of the stress field of the un- 
perturbed crack with the position of the crack edge. In 
Fourier space this contribution is given to linear order by 
HI: 5K H (k z ,x) = K°d x h+±\k z \K°^h(k z ,x) where 
v refers to the Poisson's ratio, and 8K IT (k z ,x) (resp. 

h(k z ,x)) refers to the Fourier transform of 6K n (z, x) 
(resp. h(z,x)) with respect to z. Finally, a third con- 
tribution comes from the coupling between the slope of 
the crack surface and the non singular T normal stress 
in the direction of crack propagation [2^ |2^| . This third 
contribution was shown to be negligible with respect to 
the second contribution in the thermodynamic limit |2l| . 
Finally, making the net mode II stress intensity factor 
vanish at each location z leads to: 



dh 
dx 



-J(z,x,{h}) + rj q (z,h) +r) t (z,x) 



(2) 



where the Fourier transformed elastic kernel J(fc z , x, {h}) 
is given by: 



j(k x ,x,{h}) = \k t \%-^h(k Zi x) (3) 



In other words, the morphology of the fracture surface 
h(x, z) is given by the motion of the elastic string h(z) 
that " creeps" - the x coordinate playing the role of time 
- within a random potential rj q {z,h) due to the "ther- 
mal" fluctuations rjt(z,x). The scaling properties of the 
surface h(z,x) in the steady regime are then expected 
to be described by a 2D height-height correlation given 
by Eq. ^ HH- Furthermore, if we consider an interface 
whose interaction kernel in momentum space scales as 
J(k z ,{h}) = Jo\k z \^h, the values of the critical ex- 
ponents P and z depends only on fi. In particular, 
for lon g-ra nge interaction /i = 1, one gets £ ~ 0.39 
miUll, z ~ °- 75 Ellla and = (/z ~ 0.5 in per- 
fect agreement with the values measured experimentally 
in glassy ceramics. 

The values of the critical exponents experimentally ob- 
served in homogeneous silica glass can now be discussed. 
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In this later case, the scaling properties were observed at 
length-scales smaller than the size of the process zone, i.e. 
at length-scales where the material cannot be considered 
as a linear elastic anymore. Recent AFM experiments 
[T7L fl8| and Molecular Dynamics (MD) observations [3(j 
have shown that damage spreading within this process 
zone occurs through the nucleation of nanoscale cavities 
whose nature remains controversial: They were first con- 
jectured to be similar to the ones classically observed - at 
much larger scale - during the ductile fracture of metallic 
alloys 0, . This interpretation was later questioned 
since these cavities were shown to leave no visible rem- 
nants on the post-mortem fracture surfaces [3]]]. It ap- 
pears then natural to conjecture that damage - indepen- 
dently of its precise nature - screens the elastic interac- 
tions J(z, {h}) within the process zone, making the effec- 
tive \x larger than the value y, = 1 expected for perfectly 
linear elastic materials. Rcnormalisation Groups (RG) 
methods [HH^] predict C = (2/x-l)/3, z = (5^ + 2)/9 to 
first order in e = 2/i — 1. "Arbitrary" values fi ~ 1.5 — 1.7 
would then allow us to account for the values of £ — 0.75 
and z ~ 1.2 observed in homogenous Silica glass, as well 
as for a wide range of materials 0, This has been 
confirmed through numerical simulations. Understand- 
ing how damage screening can select such an effective in- 
teraction range in crack problems provides a significant 
challenge for future investigation. 

More generally, we conjecture that both critical scal- 
ing regimes can be observed in all the heterogeneous 
materials: For length-scales smaller (resp. larger) than 
the damage zone size, ones expects Family- Viseck scal- 
ing with {C ~ 0.75, (3 ~ 0.6, z = C//3 ^ 1-2} (resp. 
{C~0.4,/?~0.5,z = C//3-0.8}). 

It should be mentioned that LEFM applied to describe 



the motion of the in-plane crack front in a disordered 
material results in a Langevin equation with non-local 
elastic kernel and quenched noise H, |28|. At the de- 
pinning transition, this approach predicts self-affine in- 
plane roughness characterised by a roughness exponent 
£ ~ 0.39 [23 and a dy namic exponent z ~ 0.75 psl |29|. 
while experiments 341] report values ( ~ 0.6 and z ~ 1. 
These experimental values are much closer to the ones 
expected in elastic line models with short range elastic 
interactionSjthat predict roughness exponents £ ~ 0.63 
and z = 1 (35j . which suggests similar damage screening 
effects as the ones invoked for out-of-plane crack rough- 
ness. 

Finally, it is worth noting that, to our knowledge, this 
study reports the first experimental observation of the 
critical roughnening predicted by the linear Langevin 
equation with non-local elastic interactions initially pro- 
posed to describe a broad variety of systems ranging 
from interfaces in disordered magnets l3o|. contact lines 
of liquid menisci on a rough substrate |37j . and inter- 
facial cracks in solids j2^. Fracture surfaces may thus 
represent an ideal experimental tool to investigate such 
pinning/depinning physics. Work in this direction is cur- 
rently under progress. 
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